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Abstract. I discuss methods and current software packages for radio searches for pulsars and 
short-duration transients. I then describe the properties of the current pulsar population and the 
status of and predictions for ongoing and future surveys. The presently observed pulsar population 
numbers around 2000 and is expected to roughly double over the next five years, with the number 
of millisecond pulsars expected to more than triple. Finally, I discuss individual objects discovered 
in the Green Bank Telescope 350-MHz Drift-Scan Survey and the Arecibo Pulsar ALFA Survey. 
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SEARCH ALGORITHMS 

Radio searches for pulsars generally rely on two types of algorithms - periodicity 
searches and single-pulse searches. The first step in either of these algorithms is 'de- 
dispersion', or correction for the frequency-dependent delays incurred as radio pulses 
travel through the interstellar medium. Filterbank data (i.e. data recorded as intensity 
vs. radio observing frequency and time) are dedispersed at a number of trial dispersion 
measures (DMs) by appropriately shifting the frequency channels in time and summing 
them to form a 'dedispersed time series'. The range of trial DMs is typically chosen 
based on a model for electron density in the Galaxy [1]. The spacing is chosen given 
some minimum detectable pulse width, with the caveat that the error due to an incorrect 
DM is not less than the dispersion smearing across an individual frequency channel. 

In a periodicity search, a Fast Fourier Transform (FFT) is applied to each dedispersed 
time series to create a power or amplitude spectrum. Multiple harmonics, typically up to 
16, are summed to increase sensitivity to narrow pulses. The time series are then folded 
at the periods of any resultant candidates to produce pulse profiles. These profiles are 
typically inspected vs. frequency and time to gauge the reality of signals. The signal- 
to-noise or x 2 values are the metrics typically used to gauge significance. In a single- 
pulse search, the dedispersed time series are analyzed by searching for any individual 
points in the series at which the radio intensity is above some threshold. The time series 
are iteratively smoothed and searched again to increase sensitivity to broadened pulses. 
Events are plotted vs. DM and time and diagnostic plots are created. Example output 
from each type of search is shown in Figure 1. For further details, see [2]. 



There are two publicly-available search packages routinely used for pulsar search- 
ing - SIGPROC 1 (with routines dedisperse and seek) and PRESTO 2 (with rou- 
tines prepdata or prepsubband and accelsearch). The de-dispersion routines 
dedisperse and prepdata use similar 'brute-force' algorithms. For large numbers 
of DMs and/or large numbers of frequency channels, prepsubband is more efficient 
as it dedisperses the data in subbands (i.e. subsets of the total bandpass) and then shifts 
and adds those subbands to form final dedispersed time series. All three routines can 
read in a list of frequency channels to ignore and have an option to clip the data to 
remove impulsive radio frequency interference (RFI). The PRESTO routine rf if ind 
will output a 'mask' file that can be read in with either prepdata or prepsubband. 
Both seek and accelsearch perform Fourier-domain periodicity searches and have 
algorithms for acceleration searches (necessary to correct for binary orbital motion). The 
accelsearch routine uses a Fourier-domain acceleration search technique [3], while 
the seek routine uses a time-domain technique that involves 'stretching' and 'squeez- 
ing' the time series [4]. Both accelsearch and seek can read in a file containing 
typical RFI frequencies to ignore in the final reported candidate list. Both search pack- 
ages also include routines for single-pulse searching. The seek routine has a single- 
pulse search option incorporated [5]. The PRESTO package has a separate routine 
named single_pulse_search . py. The algorithm used is very similar, although 
in seek adjacent pulses are added to increase sensitivity to broadened pulses whereas 
in single_pulse_search . py, time series are convolved with a boxcar function. 

CURRENT STATUS 

Figure 2 illustrates the history of pulsar discoveries. There are currently (as of January 
2011) 1880 pulsars listed in the ATNF pulsar database [6], with over 800 of these 
discovered in the Parkes Multibeam Pulsar Survey [7, 8]. Monte Carlo population 
syntheses suggest that there are roughly 120,000 detectable radio pulsars in the Galaxy 
[9], indicating that we have discovered less than 2% of all possible pulsars. The very first 
pulsars (i.e. the ~40 detected in 1968 and 1969) were discovered through their bright 
single pulses [10]. However, it was soon realized that most pulsars are detected with 
higher significance through searches for periodicity. For that reason most pulsar searches 
over the past 40 years have used solely Fourier techniques, with over 1800 of the pulsars 
listed in the ATNF pulsar database discovered through periodicity searches. Therefore, 
despite some successes (e.g. [11]), single-pulse searches were not routinely performed 
until the discovery of a new class of pulsars called Rotating Radio Transients (RRATs) in 
2006 [12]. These objects are neutron stars that have very sporadic emission that renders 
them only detectable through searches for single pulses. Giant-pulsing pulsars (e.g. [13]) 
may also be only detectable through single-pulse searches [5]. Given the small amount 
of extra computing time (typically < 10% of total) necessary for single-pulse searches, 
they are now routinely incorporated as a part of major pulsar searches. 
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FIGURE 1. Diagnostic output produced by PRESTO for periodicity (left) and single-pulse (right) 
searches. Both plots are for the same pulsar, discovered in the GBT 350-MHz Drift-Scan Survey. The 
periodicity search plot shows the folded pulse profile (top left); the folded pulse profile vs. time (bottom 
left); the folded pulse profile vs. frequency (middle top); and the x 2 of the folded pulse profile vs. DM 
(middle bottom), period (right top), period derivative (right middle), and period and period derivative 
together (right bottom). The single-pulse search plot shows the number of pulses vs. signal-to-noise (top 
left); the number of pulses vs. DM (top middle); signal-to-noise vs. DM for all pulses (top right); and all 
pulses plotted vs. DM and time, with the size of the plotted symbol corresponding to the signal-to-noise. 




1970 198D 1990 2000 2010 

Year 



FIGURE 2. The number of known radio pulsars vs. year of discovery, with a dotted line indicating 1 % 
of the total estimated detectable neutron star population [9]. 



Figure 3 illustrates the spin-down properties of the current pulsar population. Of the 
1880 pulsars currently listed in the ATNF pulsar database, 140 are in binary systems. 
Of these, ~120 have white dwarf companions and ~100 are fully recycled millisecond 
pulsars (MSPs) with periods less than 20 ms. Half of these white dwarf binary systems 
are in globular clusters. All of the Galactic systems, with the exception of the unusual 
J1903+0327 [14, 15], have circular orbits while more than 10% of the globular cluster 
systems are in eccentric (e > 0. 1) orbits. Of the remaining binary pulsars, four have main 
sequence companions and 10 have neutron star companions. The main sequence binaries 
all have very eccentric (e > 0.5) orbits while the double neutron star eccentricities range 
from 0.08 to 0.8. 

Some other classes of neutron stars are indicated on this diagram. There are over 
30 known RRATs, or objects that were discovered only through their single pulses 
[16, 17, 18, 19, 20]. Nine of those have measured period derivatives [21, 22, 23] which 
indicate higher magnetic fields than normal pulsars [21]. Magnetars, or soft gamma-ray 
repeaters and anomalous X-ray pulsars, are also shown. Three of these have been found 
to emit radio pulsations, with one of these discovered in a radio survey [24, 25, 26]. 
Finally, period derivatives have been measured for three X-ray detected isolated neutron 
stars (INSs) [27]; no radio pulsations have been detected from INSs, but this could easily 
be attributed to small beaming angles, as expected given their long periods [28]. It is 
unclear how the neutron star populations in the top-right part of the P — P diagram relate. 
One RRAT and several normal pulsars have magnetic fields comparable to magnetars. 
Two RRATs lie in a region ofP — P space closer to INSs and some lie within the bulk 
of the normal population. 

Of the 1880 pulsars listed in the ATNF pulsar database, 21 are extragalactic; all of 
these extragalactic pulsars are located in the Small and Large Magellanic Clouds. 

CURRENT SURVEYS 

Several ongoing or upcoming surveys promise to dramatically increase the known pulsar 
population and reveal interesting individual objects. We discuss the parameters of these 
surveys and the populations they probe here in order of observing frequency. Unless oth- 
erwise referenced, survey yields have been calculated using the PSRPOP software 3 [29]. 

LOFAR 140-MHz All-Sky Survey LOFAR is a wide field-of-view interferometer be- 
ing commissioned in the Netherlands, with over 30 antennas currently operational, 
including one of the UK and one in Sweden [30]. At a very low center frequency 
of 140 MHz, with 26 MHz of bandwidth and 64-min integrations, an all-sky sur- 
vey would discover roughly 900 pulsars [31], including nearly all of the detectable 
pulsars within ~1 kpc. LOFAR will also be a powerful probe of pulsars in globular 
clusters and local group galaxies. 

GBT 350-MHz Drift-Scan Survey Over the (Northern) summer and fall of 2007, 
10,300 deg 2 between declinations of -21 and +26 were covered in 'drift-scan' 
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FIGURE 3. Period vs. period derivative for all known pulsars. Pulsars are plotted as dots, with those 
discovered in the PALFA/GBT 350-MHz Drift-Scan Survey plotted as diamonds/triangles. Binary pulsars 
are marked with circles, with the ellipticity of the circle corresponding to the ellipticity of the binary orbit. 
RRATs are denoted by open stars, SGRs by plus signs, and INSs by crosses. Lines of constant inferred 
surface dipole magnetic field strength and characteristic age are shown. 



mode during the Green Bank Telescope (GBT) track repair, with an effective inte- 
gration time of 150 s. A center frequency of 350 MHz, 50 MHz bandwidth, 2048 
frequency channels, and 81.92 /is sampling time were used. A total of 26 pulsars, 
including five MSPs and one RRAT, have been discovered, with processing roughly 
80% complete 4 . A portion of the data (20 Terabytes, or 1500 deg 2 ) has been set 
aside for analysis by students in the Pulsar Search Collaboratory (PSC) project 
[32]. The students have discovered two pulsars and several more are expected. 
GBT 350-MHz Northern Celestial Cap Survey The primary motivation of this sur- 
vey is to find MSPs to include in the pulsar timing array for gravitational wave de- 
tection. This GBNCC survey has recently begun and will cover declinations greater 
than 38° at 350 MHz with two-min integrations using 100 MHz bandwidth, 4096 
frequency channels, and 81.92 /is sampling time. It is expected to discover over 
200 normal pulsars and several tens of MSPs. 



http : / / www . as . wvu . edu/ -pulsar/ GBTdrif t350/ 



GMRT 330- and 610-MHz Surveys Three pulsars have been discovered in a prelimi- 
nary survey with the Giant Metrewave Radio Telescope in India covering \b\ < 1° 
and 45° < I < 135° [33] and using a frequency of 610 MHz, with 16-MHz band- 
width, 256-/1S sampling time, and 35-min integrations. Up to 100 pulsars (including 
~5 MSPs) could be discovered in a planned 350-MHz survey with 32-MHz band- 
width 1° < \b\ < 10°. The GMRT has the ability to image survey fields-of-view 
concurrently with the search observations, facilitating pulsar localization. 

Arecibo 1.4-GHz Pulsar ALFA Survey This ongoing seven-beam survey covers \b\< 
5° and 40° < / < 75°, 170° < / < 210° with 300 MHz of bandwidth, 256 frequency 
channels, 64-/is sampling time, and 265-s integrations [34]. So far, 56 pulsars, 
including eight MSPs and six RRATs, have been found 5 . Roughly 300 pulsars, 
including about 50 MSPs, are expected to be discovered. 

Parkes 1.4-GHz High Time Resolution Universe Survey This survey will cover the 
entire sky visible from Parkes with 340-MHz bandwidth, 1024 frequency channels, 
and 64-jUs sampling time. Integrations of 4300, 540, and 270 s will be used at 
\b\ < 3.5°, 3.5 < \b\ < 15°, and high latitudes, respectively. The mid-latitude part 
of the survey has begun and has discovered 27 pulsars, including five MSPs [35]. 
The full survey is expected to discover ~400 pulsars, including 75 MSPs. 

Effelsberg 1.4-GHz Multibeam Survey This survey will cover the same latitude 
ranges as Parkes with 1500 s, 3 min, and 1.5 min integrations with a seven-beam 
receiver identical to that used for PALFA. The survey will use a bandwidth of 240 
MHz, 1024 frequency channels, and 54-/is sampling time. Roughly 600 pulsars 
are expected to be discovered, with ~100 likely to be MSPs [36]. 

Parkes 6-GHz Galactic Center Survey This recently completed survey covered \b\ < 
0.25° and —60° < / < 30° using a seven-beam receiver with a 576-MHz bandwidth, 
192 frequency channels, and 125-/1S sampling time and integration times of 1055 s. 
The survey discovered three pulsars, all with DMs greater than 900 pc cm~ 3 [37]. 

A total of 4500 pulsars are predicted to be discovered by the above surveys combined. 
Even accepting that half of these pulsars will be discovered in multiple surveys, we 
still expect the pulsar population to roughly double over the next five years. This will 
enable precise population analyses, studies of pulsar emission phenomenology, and a 
clearer understanding of how radio pulsars are related to other classes of neutron stars. 
Given the numbers of RRATs detected in current surveys, we expect the number of these 
objects to increase dramatically, to at least several hundred objects. This will allow us 
to perform population analyses to clarify how they are related to normal pulsars. The 
Galactic MSP population, which currently numbers around 100 (including the Fermi 
discoveries, e.g. [38]), will more likely triple over the next five years. This will nearly 
double the sensitivity of the pulsar timing array for gravitational wave detection [39]. 
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GBT 350-MHZ DRIFT-SCAN SURVEY 



The discoveries from this survey are marked in Figure 3 and include J1023+0038, the 
first MSP showing evidence for recent accretion [40] and a unique probe of the MSP re- 
cycling process. Another highlight is the 'black-widow' pulsar J2256-10, a 2-ms eclips- 
ing pulsar with a minimum companion mass of only 0.03 M [41]. We have recently 
detected gamma-ray pulsations with Fermi from this source. A third notable discovery 
is J2222-01, a 33-ms partially recycled pulsar with a massive (i.e. minimum mass of 
1.1 M ) white dwarf companion [23]. This object seems to belong to the rare class of 
intermediate mass binary pulsars [43]. It has a DM of only 3.27 pc cm~ 3 , yielding an 
inferred distance [1] of only 300 pc. This is the second closest pulsar binary and will be 
an excellent target for high-precision astrometry and multi-wavelength companion char- 
acterization. A final interesting object is J0348+04, which is only partially recycled (i.e. 
period of 39 ms) yet has a very low mass companion (i.e. minimum mass of 0.08 M ) 
[42]. This object is similar to pulsar J1744-3922, discovered with Parkes [44]; together 
they may represent a new class of neutron star binary system. 

PULSAR ALFA SURVEY 

The discoveries from this survey are also marked in Figure 3 and include the intriguing 
J1903+0327, a fully recycled (2.1 ms) MSP in an eccentric (e = 0.44) orbit with a 
solar mass companion [14]. The companion is a main-sequence star, suggesting that 
the pulsar was originally a member of a triple system in which the donor star was either 
ablated or ejected [15]. The pulsar has a very high and unusually well-determined mass 
of 1.667±0.021 M . The survey has also discovered a new double neutron star system, 
J 1906+0746 [45], in which the younger, non-recycled pulsar is the detectable source. 
This is the youngest binary pulsar system known. Seven pulsars have been discovered 
through the single-pulse search, including one mysterious source (J 1928+ 15) which 
showed three consecutive detectable pulses but none since, despite multiple follow-up 
attempts [19]. Some of the survey data are being analyzed through the Einstein® Home 
distributed computing project 6 [46]. This has so far resulted in the discovery of two new 
pulsars. Both of these pulsars are recycled and one may have a neutron star companion. 

ACKNOWLEDGMENTS 

I gratefully acknowledge my collaborators on the GBT 350-MHz Drift-Scan and PALFA 
survey teams, Duncan Lorimer for prediction of survey yields using the PSRPOP soft- 
ware, and the ASTRONS2010 organizers for the conference invitation and the travel 
support. This work has been supported by WVEPSCOR, the Sloan Foundation, the Re- 
search Corporation, and the NSF 



http: / /einstein . phys . uwm . edu 



REFERENCES 



1. J. M. Cordes, and T. J. W. Lazio, astro-ph/0207156 (2002). 

2. D. R. Lorimer, and M. Kramer, Handbook of pulsar astronomy, Cambridge observing handbooks for 
research astronomers, Vol. 4. Cambridge, UK, Cambridge University Press (2005). 

3. S. M. Ransom, S. S. Eikenberry, and J. Middleditch, Astron. J. 124, 1788 (2002). 

4. F. Camilo et al, Astrophys. J. 535, 975 (2000). 

5. J. M. Cordes, and M. A. McLaughlin, Astrophys. J. 596, 1 142 (2003). 

6. R. N. Manchester, G. B. Hobbs, A. Teoh, and M. Hobbs, Astron. J. 129, 1993 (2005). 

7. R. N. Manchester et al, Mon. Not. R. Astron. Soc. 328, 175 (2001). 

8. R. P. Eatough et al, Mon. Not. R. Astron. Soc. 407, 2443 (2010). 

9. C.-A. Faucher-Giguere, and V. M. Kaspi, Astrophys. J. 643, 332 (2006). 

10. A. Hewish et al, Nature, 217, 709 (1968). 

11. D. J. Nice, Astrophys. J. 513, 927 (1999). 

12. M. A. McLaughlin et al, Nature 439, 817 (2006). 

13. J. M. Cordes et al, Astrophys. J. 612, 375 (2004). 

14. D. J. Champion et al, Science 320, 1309 (2008). 

15. P. C. C. Freire et al, arXiv: 101 1.5809 (2010). 

16. R. P. Eatough, E. F. Keane, and A. G. Lyne, Mon. Not. R. Astron. Soc. 395, 410 (2009). 

17. E. F. Keane et al, Mon. Not. R. Astron. Soc. 401, 1057 (2010). 

18. S. Burke-Spolaor, and M. Bailes, Mon. Not. R. Astron. Soc. 402, 855 (2010). 

19. J. S. Deneva et al, Astrophys. J., 703, 2259 (2009). 

20. J. W. T. Hessels et al, "The GBT350 Survey of the Northern Galactic Plane for Radio Pulsars and 
Transients," in 40 Years of Pulsars: Millisecond Pulsars, Magnetars and More, edited by C. Bassa, 
Z. Wang, A. Cumming, and V. M. Kaspi, vol. 983 of American Institute of Physics Conference Series, 
613 (2008). 

21. M. A. McLaughlin et al, Mon. Not. R. Astron. Soc. 400, 1431 (2009). 

22. A. G. Lyne et al, Mon. Not. R. Astron. Soc. 400, 1439 (2009). 

23. J. Boyles et al, in preparation (201 1). 

24. F. Camilo et al, Nature 442, 892 (2006). 

25. F. Camilo, S. M. Ransom, J. P. Halpern, and J. Reynolds, Astrophys. J. 666, L93 (2007). 

26. K. Lazaridis et al, Mon. Not. R. Astron. Soc. 390, 839 (2008). 

27. D. L. Kaplan, and M. H. van Kerkwijk, Astrophys. J. 692, L62 (2009). 

28. V. I. Kondratiev et al, Astrophys. J. 702, 692 (2009). 

29. J. P. Ridley, and D. R. Lorimer, Mon. Not. R. Astron. Soc. 404, 1081 (2010). 

30. J. Hessels et al, arXiv: 1009.1758 (2010). 

31. J. van Leeuwen, and B. W. Stappers, Astro. & Astroph. 509, A7 (2010). 

32. R. Rosen et al, arXiv: 1005. 1060 (2010). 

33. B. C. Joshi et al, Mon. Not. R. Astron. Soc. 398, 943 (2009). 

34. J. M. Cordes et al, Astrophys. J. 637, 446 (2006). 

35. M. J. Keith et al, Mon. Not. R. Astron. Soc. 409, 619 (2010). 

36. D. R. Lorimer et al, Mon. Not. R. Astron. Soc. 372, 111 (2006). 

37. S. D. Bates et al, Mon. Not. R. Astron. Soc. 1820 (2010). 

38. S. M. Ransom et al, Astrophys. J. 727, L16, (201 1). 

39. G. Hobbs et al, Classical and Quantum Gravity 27, 084013 (2010). 

40. A. M. Archibald et al, Science 324, 141 1 (2009). 

41. I. H. Stairs et al, in preparation (201 1). 

42. R. Lynch et al, in preparation (201 1). 

43. R. D. Ferdman et al, Astrophys. J. 711, 764 (2010). 

44. R. P. Breton et al, Astrophys. J. 661, 1073 (2007). 

45. D. R. Lorimer et al, Astrophys. J. 640, 428 (2006). 

46. B. Knispel et al, Science 329, 1305 (2010). 



